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A dedicated cosmic muon Monte-Carlo event generator CMSCGEN has been developed for the 
CMS experiment. The simulation relies on parameterisations of the muon energy and the in- 
cidence angle, based on measured and simulated data of the cosmic muon flux. The geometry 
and material density of the CMS infrastructure underground and surrounding geological layers 
are also taken into account. The event generator is integrated into the CMS detector simulation 
chain of the existing software framework. Cosmic muons can be generated on earth's surface as 
well as for the detector located 90 m underground. Many million cosmic muon events have been 
generated and compared to measured data, taken with the CMS detector at its nominal magnetic 
field of 3.8 T. 
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1. Muon flux parameterisation and normalisation 

The cosmic muon generator of the CMS experiment [|I|] is based on a parameterised differential 
muon flux obtained by means of the air shower program CORSIKA ^ which in turn has been 
operated making use of the EPOS interaction model for high energy interactions and the GHEISHA 
interaction model for low energy hadronic interactions. The CORSIKA simulation results have 
been fitted with polynomials [Q] to describe the differential flux at earth's surface as a function 
d J dcc^e d(t) °^ ^^^ azimuth angle (j), the incidence angle 6 and the muon momentum p. The flux 
is normalised to vertical muons from measured data with an energy of 100 GeV. The maximally 



allowed phase space in azimuth angle is G [0, 27r), whereas azimuthal isotropy has been assumed. 
The phase space for incidence angles is given by 0° < < 84°. For incidence angles of 6 > 75° the 
parameterisation is extrapolated. In momentum the phase space is restricted to 3 < p^ < 3000 GeV. 
The lower limit is driven by the existence of new physical processes which set in at such low 
momenta which makes the extrapolation invalid. The upper limit of muon momenta is driven by 
strongly decreased flux in this regime. The muon momentum is approximated by a polynomial of 
the expression 

L = logio(WGeV) (1.1) 

which is slowly varying in muon momentum. The muon momentum spectrum is then fitted by a 
polynomial of the form 

s{L)=ao+aiL + ...+a(,L^. (1.2) 

The momentum dependent zenith angle is taken into account by a polynomial of the cosine of the 
incident angle 

z(cos d,L)= bo{L) + bi (L) cos d + b2{L) cos^ d. (1.3) 

Finally the differential flux is given by 

d(i> 1 1 

Cnorm-^-^(L)-z(cOS0,L)-— . (1.4) 



dp dcos 6 d^ p^ '271 

2. Muon energy loss in material 

A random vertex on a disk at the surface is assigned to a generated single muon. The size of 
the disk which is centred around the vertical symmetry axis of the detector is determined by the 
chosen maximal incidence angle such, that a muon whose extrapolated direction does hit the CMS 
detector does necessarily have a vertex inside the disk. In the case the muon direction is pointing to 
the CMS detector the muon is propagated through the different material densities of the geological 
environment between the surface and the CMS detector, located at 90 m underground to obtain an 
integrated amount of water equivalents. This amount determines the energy loss of a given muon 
as a parameterised function [^ of the energy of the incident muon at surface. If the direction 
is not pointing to the CMS detector a new random vertex at the surface is chosen. A material 
map describes the diverse materials from earth's surface to the CMS detector. In concrete, the 
fundament of the hall at surface, the three vertical access shafts, a movable plug of the main shaft 
as well as the collision and service caverns including the adjacent parts of the LHC tunnel are taken 
into account. Two different average densities are assigned to the geological layers surrounding the 
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Figure 1: Left plot: Entrance vertex at the surface of the generated muons which arrive at the CMS 
detector (indicated by the dashed rectangle) 90 m underground. The three vertical shafts are indicated by 
black circles. There is clearly a correlation between muon intensity and shafts visible. Right plot: Energy 
spectrum of muons in data (black points) measured with the CMS detector In comparison the GEANT 
simulation (green/blue histogram) is shown. The simulation has been normalised to the number of events of 
the data. 

CMS infrastructure. The upper half consists of sand, clay, gravels and water while the lower half 
consists of rock. In Fig. [T] the vertex position of the muons at the surface which reach the CMS 
detector are shown. Lower energetic muons are responsible for the enhanced intensities at the 
vertical access shafts. After the generated muons have been propagated to the CMS detector a full 
GEANT [^ simulation of the detector processes the muons further. 

3. Simulation and data comparison 

In the cosmic data taking period 2008 three hundred million events have been recorded with a 
solenoid magnetic field of B = 3.8 T. Global cosmic muons, which are reconstructed in the muon 
chambers and the central tracking system are chosen. A mixture of four different triggers is used 
in the simulation. The momentum of the reconstructed global muons is required to exceed a mo- 
mentum threshold of /? = 10 GeV. The incidence angle is restricted to the interval < < 60° to 
ensure that the parametrisation is applied far away from its regime of extrapolation. Fig. |] right 
plot shows the energy spectrum of the cosmic muon data taken in 2008 in black points. The sim- 
ulation is superposed as a green/blue shaded histogram and normalised to the number of entries 
in data. The decrease of the data distribution by the power of -2.7 is very well described by the 
simulation. 
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